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ABSTRACT 

The nature of ultraluminous X-ray sources (ULXs) is presently unknown. A possible explanation is 
that they are accreting intermediate mass black holes (IBHs) that are fed by Roche lobe overflow from 
a tidally captured stellar companion. We show that a star can circularize around an IBH without 
being destroyed by tidal heating (in contrast to the case of M, > 10 6 Af Q massive black holes in galactic 
centers, where survival is unlikely). We find that the capture and circularization rate is ~ 5xl0~ 8 yr , 
almost independently of the cluster's relaxation time. We follow the luminosity evolution of the 
binary system during the main sequence Roche lobe overflow phase and show it can maintain ULX- 
like luminosities for > 10 7 yr. In particular, we show that the ULX in the young cluster MGG-11 in 
star-burst galaxy M82, which possibly harbors an IBH, is well explained by this mechanism, and we 
predict that > 10% of similar clusters with IBHs have a tidally captured circularized star. The cluster 
can evaporate on a time-scale shorter than the lifetime of the binary. This raises the possibility of a 
ULX that outlives its host cluster, or even lights up only after the cluster has evaporated, in agreement 
with observations of host-less ULXs. 

Subject headings: black hole physics — stellar dynamics — galaxies: star clusters — X-ray binaries 



1. INTRODUCTION 

Black holes (BHs) have deep potential wells and can 
transform gravitational energy very efficiently to radi- 
ation. The energetic central engines of quasars are 
thought to host massive black holes (MBHs) of M. > 
10 6 M Q . It is natural to extrapolate this idea and invoke 
an intermediate mass black hole (IBH; 10 2 < M./M Q < 
10 5 ) to explain ultraluminous X-ray sources (ULXs), 
which are considerably brighter than a stellar mass ob- 
ject radiatin g at its Eddington luminosity. For example, 
Kaaret et al. If20f)lt) have suggested that an IBH powers 
the ULX in MGG-11 in star-burst galaxy M82. 

The origin of the gas that fuels the X-ray source is un- 
clear, since almost all the gas in young clusters is rapidly 
blown away by the strong winds of massive stars. One 
possibility for providing the gas is the tidal disruption of 
a main sequence (MS) star of mass M* and radius i?* 
with periapse r p < r t , where r t = (M./M*) 1 / 3 /?* is the 
tidal radius. However, direct disrupt ions le ad to a short 
flare (t fl „ rK < yr; Ree s! 19881: Ulme J 19991 Ayal, Livio, 
& Pirar 2000), which is incompatible with the ~ 20 yr 
observation period of the X-ray source. In this Letter we 
investigate a more gradual process for feeding the IBH, 
namely the tidal capture of a MS star and the subsequent 
Roche lobe overflow. 

2. TIDAL CAPTURE RATE 

A BH in a cluster with velocity dispersion a dominates 
the potential within its radius of influence r a = GM m /a 2 ; 
inside r a orbits are approximately Keplerian, and stars 
are distributed according to a power law n* oc r~ Q , with 
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a w 3/2 (Bahcall & Wol fl 19761 Baumgardt, Makino, 
& Portegies Zwar f2003l) . The cusp is truncated inside 



some radius 



e.g. 



(M,/M*)R* where the rate 



of destructive colli sions e xceeds the two-body relaxation 
rate (Frank & Ree i 197fij) . 

Stars can reach an orbit with periapse of order of the 
tidal radius by angular momentum diffusion. When the 
star passes at r p , an energy AE t (r p ) is invested in raising 
tides, causing the star to spiral in (r p < 3r t is typically 
required for an appreciable effect). The evolution of a 
tidally heated star is not well-understood. Two extreme 
models of "squeezars" (stars that are continually pow- 
ered b y tidal squeezing) were studied by Alexander & 
Morris f?200l . "Hot squeezars" are heated only in their 
outer layers and radiate their excess energy efficiently; 
they hardly expand. "Cold squeezars" dissipate the tidal 
energy in their bulk and puff up to giant size. 

Our analysis is based on the following assumptions. 

(1) The stars are "hot squeezars" (in we discuss 
some consequences of relaxing this assumption). 

(2) As long as the eccentricity e is high, 



1 



,/a < £ e ~ 0.1, 



(1) 



where a is the orbital semi-major axis, the stellar struc- 
ture is not significantly affected by the tidal heating, and 
the tidal energy dissipated per orbit, 



AE t {b) = 



GUI T(b) 
R* b e 



(2) 



is constant (Alexander & Morri f2003fl : here b = r p /r t , 
and T(b) is the tidal coupling coefficient, which depends 
on the stellar structure and is a s trongl y decreasing func- 
tion of b (e.g. Press & Teukols H 19771) . When the orbit 
decays to the point where 1 — e > £ e , the tidal heating 
drops off until e ventual ly the star circularizes at a > r t 
and AE t = (Hu£TM3). 

(3) The star can survive as long as its tidal luminosity 
does not exceed, to within order unity, its Eddington 
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luminosity L E = 1.3 x 10 38 ergs -1 M*/M , 

AE t /P < £ L L E , (3) 

where P is the orbital period and 

The tidal heating rate is highest when P is shortest, 
just before tidal heating shuts off when a = 6r t /£ e (Eq. 

Therefore, the Eddington luminosity limit (Eq. |3J) 
corresponds to a minimal periapse bminTt that a star can 
have and still circularize without being disrupted, which 
is given implicitly by 



2- 



in 



3/2 



(4) 



When r p < b min r t , the star is evaporated by its own 
tidally powered luminosity during in-spiral. 

Stars within the "loss-cone" , a region in phase space 
where stars have periapse s maller than rt (Frank & Rees 
I 197ft Lightman & ShaDirc Tl977t) . are disrupted by the 
BH. Two-body scattering sustains a flow of stars in an- 
gular momentum space towards the loss-cone. During 
in-spiral, two-body interactions change the periapse of 
the star. The time t p over which the periapse of a star is 
changed by order unity due t o many small angle deflec- 
tions is (Alexander & Hopmar 2003) 



t p (b, a) = —t r , 
a 



(5) 



where t r is the relaxation time. Note that t r does not 
depend on the distance from the BH for a — 3/2. 

The in-spiral time is the time it takes until the semi- 
major axis of the star becomes for mally z ero; for a hot 
squeezar it is (Alexander & Hopmar. 20 03() 



t (b,a) 



2-irMWGM.a 
AE t (b) 



(6) 



If deflections increase the periapse, the dissipation be- 
comes much less efficient, while if the periapse decreases, 
the star may cross r t and be disrupted. Either way 
it fails to circularize. Circularization can happen only 
if the in-spiral time to is shorter than the time-scale 
for deflections, t p . The widest orbit a c (b) from which 
a star can still spiral in for periapse br t i s give n by 
to{b,a c ) = 3t p (b 7 a c ) (Alexander & Hopman 2003, Eq. 
[11] for a = 3/2). It then follows from Eq. © that the 
maximal distance a max from which a star can originate 
to reach the tidal radius without being destroyed is 



3A£ , t (6 m in)6 m i n rttr 



2nM±y/GM. 



2/3 



(7) 



Within r in the cusp flattens and relaxation is inefficient, 
so there are hardly any stars on eccentric orbits. Since 
ri n grows more rapidly with M, than a max , there exists 
a maximal BH mass M max , such that for M. > M max , 
amax < Hn, and tidal capture is strongly suppressed. Fig. 
CJ shows a max and r- ln as a function of M # ; for the cal- 
culation of t r we assumed that the M, — a relation 



M. = 1.3 x 10 M G 



( Z V 

1 200 kms- 1 / 



(8) 



,200 kms 

(Ferrarese & Merrit t HEXXl Gebhardt et al. 1^00(1 
Tremaine et al J 2002j) can be extend ed to IBHs (see e.g., 
Portegies Zwart & McMillan ! 20 02). Circularization is 
only possible for M, < M max « 1O 5 M . 




log(M/Msun) 

Fig. 1. — Dependence of a max (solid), r In (dashed), and rt 
(dashed-dotted) on the mass of the BH, for a IOMq star. Cir- 
cularization is only possible provided that a max > r; n . 



The rate T at which stars diffuse into orbits that al- 
low s uccessf ul circularization is given by (Eq. [9] Syer & 
Ulme fl999l) 



r = 



(flmax/ra)^^, 



Mn(2- v /a max /i>min»''t) 



(dmax > r in ), (9) 



where the logarithmic term expresses the depletion of 
the stellar density near the loss-cone; N a is the number 
of stars within the radius of influence. The rate is essen- 
tially independent on M+ for a fixed stellar mass within 
r a , and it decreases only logarithmically with t r (cf . Eqs 
|3 15] ) : a larger t r increases the volume of stars that con- 
tributes to T, but decreases the rate at which stars enter 
the loss-cone. The rate does not depend very sensitively 
on our assumptions: roughly, T oc ^l/^c- 

3. ROCHE LOBE OVERFLOW ON THE MAIN SEQUENCE 

Orbital angular momentum conservation implies that 
the circularization radius is a c i rc = 2b m { n r t . Efficient in- 
spiral and successful circularization require & m i n ^2 — 2.5, 
so that a c ; rc ~ (4 — 5)r t . The onset of mass transfer 
through the Roche lobe occurs when the distance be- 
tween the IBH and the star is a c i rc ~_2rv (assuming 
M* = 10M Q , M*/M. - 0.01, EggletorJ 1983i). This is 
roughly a factor of two smaller than the typical value of 
Ocirc- However, as it evolves on the MS, a 1OM star 
expands by a factor of ~ 2.7 by the time it reaches the 
terminal age MS (TAMS). This implies that Roche lobe 
overflow (RLOF) occurs at some point on the MS, and 
continues for a time t x , which is shorter than the MS 
lifetime £ms (we assume here that the star expands sig- 
nificantly only after it has circularized). In the following 
analysis we assume for simplicity that RLOF holds over 
the entire MS. This does not affect our conclusions sig- 
nificantly since the observationally relevant phase of high 
X-ray luminosity occurs toward the TAMS. For massive 
MS stars RLOF is preceded by a less luminous phase 
resulting from the accretion of strong stellar winds. 

Mass transfer from a MS star to an IBH is driven by the 
thermal expansion of the donor and the loss of angular 
momentum from the binary system. Mass transfer then 
implies that the donor fills its Roche lobe (i?* = Rm) 
and continues to do so (R+ = i?Ri). We assume that as 
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log(t/Myr) 

Fig. 2. — Mass evolution of the donor star, assuming Roche-lobe 
contact at the zero-age MS to a 1400A/q IBH. The solid, dashed 
and dotted curve are for a 15Mq, IOMq, and 5Mq donor. 



long as the Eddington limit is not exceeded, all the mass 
lost from the donor via the Roche lobe is accreted by the 
IBH (M. = —Mi,). Otherwise, the mass in excess of the 
Eddington limit is lost from the binary system. 

The expansion of the donor on the MS is ca lculated us- 
ing fits from Eggleton, Tout & Fitchett (1989) to detailed 
stellar evolution calculations. We assume that the evo- 
lution of the donor was not affected by mass loss. Varia- 
tions in the Roche radius of the donor can be computed 
from the redistribution of mass and angular momentum 
in the binary system. The radius of the Roche -lobe i s 
estimated with the fitting formula from Eggleton) (1983). 

We stop following the binary evolution at the TAMS; 
the simple model for calculating the amount of mass 
transfer may be inappropriate for the post-MS evolu- 
tion of the donor, as the star then rapidly expands. 
However, at the end of the MS the donor still has a 
considerable envelope and the star ascends the giant 
branch. The post-MS evolution is likely to result in a 
short (£pms < 0.1 £ms) phase in which the luminosity 
increases by more than an order of magnitude. 

In Fig. we plot the mass of the donor as a function 
of time. It is assumed here that there is Roche lobe 
contact directly after circularization. As discussed, this 
actually only happens when the star has evolved towards 
a later stage of the MS phase, which is when the mass 
loss in the plot starts to drop more rapidly. 

We estimate the X-ray luminosity during mass transfer 
with t he model discussed by Kording, Falcke, & Markoff 
(2002). They argue that the X-ray luminosity is gener- 
ated by an accretion disk. The binary is in the hard state 
if M > M cr it, in which case L x = eMc 2 . At lower accre- 
tion rates L x = eMc 2 M / M cr n, in which case the X-ray 
source becomes transient (i.e. short outbursts, separated 
by long states of quiescence, Kalogera et al. 2003). For 
Mfrit we adopt the equation derived by Dubus et al. ( 
1999, see Eq. 32) and assume e = 0.1. These choices 
are co mparab le to M crrt ~ 10 _7 M Q yr _1 of Kording 
et al. L(2QQ2l) . The resulting X-ray luminosity is pre- 
sented in Fig. ©. Note that lower mass donor binaries 
(M* ,$ 5M Q ) live longer, are less luminous and tend to 
show transient behavior, where high mass donor binaries 
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Fig. 3. — X-ray luminosity for a 1400Mq accreting IBH as a 
function of time. Line styles as in Fig.|3 To the right side of 
the figure we added M in logarithmic units of Mgyr" 1 for the 
regime where M > M cr it . The radius of the star grows significantly 
towards the TAMS (where the lines for 15Mq and 10 Mq donors 
rise). It is probably only near that point at which RLOF actually 
starts, so that the luminosity is higher than would be estimated 
from RLOF of a zero age MS star. The 5 Mq donor does not show 
this rise in L x as M < M cr ; t after ~ 10 Myr. This also causes the 
drop of L x for the IOMq donor star at the end of its evolution. 
When M < M cr it, the source becomes transient. 



are more luminous, shorter lived, steady sources. 

4. CLUSTER MGG-11 

We apply our analysis to the young dense star clus- 
ter MGG-11 in the irregular galaxy M82, at a distance 
of ~ 4 Mpc. This cluster contains the variable X-ray 
source M82-X7 with L x = (0. 8-160) x 10 39 ergs" 1 (Wat- 
son, Stange r fc Gr iffith ! 19841 Matsumoto & Tsur A 199ft 
Kaaret et al l 2001|) . The velocity dispersion in the cluster 
is accurately meas ured, a = 11.4±0.8 kms" 1 (McCrady, 
Gilbert, & Grahar J2003^ . We assume that this cluster 
contains an IBH which is also the engine for the X-ray 
source. If this IBH obeys the M. — a relation (JSJl, its 
mass is M, — 1.4 x 10 3 A/ Q , which is consiste nt wit h 
the recent calculations of Portegies Zwart et al. [(2004), 
who show that an IBH could have formed dynamically 
in MGG-11 by a runaway merger of MS stars. Within 
its radius of influence r a = 0.05 pc, the number of stars 
is N a = 2M./M+ (MerritOooS) and t r ~ 10 5 - 10 6 yi\ 

The age of the cluster is t c i = (7 — 12) Myr, corre- 
sponding to a turn-off mass of 17 — 25M Q (Eggleton et 
al l 19891) . The mean stellar mass of the cluster at birth is 
(A/*) = 3 Mq, but as a result of mass-segregation the av- 
erage mass within r a is much higher. The dir ect N-body 
calculations of Portegies Zwart et al. 1(20041) show that 
at an age of 7 Myr, the mean mass of the single stars 
in the core of MGG-11 is (M*) = 8 ± 3M . For sim- 
plicity we assume within r a a single mass population of 
stars with M + = 10 Mm, and radius R+ — 5.4i? Q (Gorda 
& SvechnikoA|j998). The results of Fig. show that 
a MS donor of mass M+ > 10 Mq can account for the 
luminosity of the ULX in MGG-11. 

We assume (CcCl) = (0.1,0.5), and take the numeri- 
cal values for the fun ction T (b) for parabolic orbits from 
Alexander & Kumar I? 20011) . With these parameters we 
find a capture rate of T = 5 x 10 _8 yr . This implies 
that a fraction T t c it x /tMS — (30 — 50)% of clusters har- 
boring an IBH has formed a tidal binary and may be 
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observed during RLOF in the MS phase. A fraction 
rtcZ^PMs/iivis ~ 4% of clusters with an IBH should be 
observed during the more luminous post-MS phase. 

5. SUMMARY AND DISCUSSION 

MS stars can spiral into an IBH as a result of tidal cap- 
ture and circularize close to the tidal radius. This pro- 
cess is unique to IBHs, since stars cannot survive tidal 
in-spiral around a MBH in a galactic center. After cir- 
cularization, the star expands on the MS until high lu- 
minosity RLOF accretion starts toward the end of the 
MS phase. We analyzed RLOF during the MS phase in 
some detail and calculated the X-ray luminosity. Post- 
MS RLOF is harder to model, but the resulting lumi- 
nosity is expected to be at least an order of magnitude 
brighter and about an order of magnitude shorter in du- 
ration. The X-ray luminosity is consistent with observed 
ULXs, such as MGG-11. 

MGG-11 is the only cluster out of hundreds in M82 
with a ULX. Possibly other clusters we re not sufficiently 
dense to fo rm IB Hs (Matsushita et al.£.2000; Portegies 
Zwart et al l 2004|) . If a fraction /. of the N c i clusters in 
M82 harbors an IBH, the number of ULXs is estimated 
by iV x = f,Tt c it x /tMS- Thus /. has to be of the order of 
a few percent in order to account for one ULX in M82. 

In order to circularize, a star has to dissipate ~ 
(M./M*) 2 ' 3 times its binding energy. If a certain frac- 
tion 5 of the energy is invested in bulk heating (for "hot 
squeezars" S = as assumed so far, for "cold squeezars" 
(5=1), the star expands. An X-ray binary can form only 
if S < (M+/M.) 2/3 . Nevertheless, a shorter lived ULX- 
phase is still possible even if <5 > (M*/M.) 2 / 3 . When the 
star expands to a radius > bR+ it is gradually peeled ev- 
ery periapse passage and feeds the IBH for a period much 
longer than ifl arc . However, the process is limited by the 



two-body deflection time-scale t p , which is < 10 3 yr for 
a 10 3 Af Q IBH. This translates to a detection probability 
of only rtp~5xl0 -5 , and so it is unlikely that the ULX 
in MGG-11 originates in this type of process. For simi- 
lar reasons, it is improbable to observe a very luminous 
tidally heated star (squeezar) during the final stages of 
its in-spiral into an IBH in a stellar cluster (this may 
be possible for squeezars near the MBH in the Galactic 
Center, where the in-spiral time is long er and the capture 
rate higher, Alexander & Morris 20Q3). 

The lifetime of the host cluster is limited by the galac- 
tic tidal fie ld, and can be as short as 100 Myr (Portegies 
Zwart et al l 20011) . This is much shorter than the RLOF 
phase of a low mass donor (e.g. ~ Gyr for a 2M star). 
Thus, the X-ray life-time of a low- luminosity binary can 
be much longer than the life-time of the cluster. Qui- 
escent orphaned IBHs can suddenly light up when their 
companion ascends the giant branch and starts to trans- 
fer mass to the IBH. Our scenario predicts the existence 
of host-less ULXs, which are more likely to be transient 
and less luminous. Their exact fraction in the ULX pop- 
ulation cannot be reliably estimated at this time. How- 
ever, it is interesting to note that 3 — 10 out of 14 of 
the ULXs in the Antennae Galaxies are coincident with 
a stel lar cluster, while the others are not (Zezas et al. 

cm. 
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